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We study the binding energy, intertube distance and electronic structure of bundles consisting 
of single walled carbon nanotubes of the same chirality. We model various nanotube structures 
(chiralities) and orientations with van der Waals density functional theory. The orientation of the 
tubes in the bundle strongly influences the properties of the bundles if the chirality of the tubes 
shares symmetry with the trigonal bundle structure, meaning chiralities which have a 60°-rotational 
symmetry (C6-axis), e.g. (12,0) bundles. The bundle structure breaks the symmetry depending 
on the arrangement of the neighboring tubes. Pseudogaps open in the electronic density of states 
and intertube distances (±5-10%) vary in dependence of the relative orientation of the tubes in 
the bundle. Bundles of C6-axis armchair tubes have metallic configurations. A 15° rotation off 
the high symmetry configuration (AA-stacked) of a (6, 6) bundle shows metallic behavior and has 
higher binding energy than the high symmetry configuration. We find binding energies between 
19 meV/atom and 35 meV/atom, depending on the chirality of the tubes. The intertube distances 
are between 3.2 A and 3.4 A but independent of orientation for non C6-axis tubes. 

PACS numbers: 61.46.-w, 61.48.De, 71.15.Nc, 71.20.-b, 73.22.Gk 



I. INTRODUCTION 

Carbon is one of todays most exciting materials with 
properties originating from the structure and symme- 
try of its allotropesii~— The properties of carbon nan- 
otubes depend on their one dimensionality and the ex- 
act arrangement of the carbon atoms on their surface, 
called chirality 2,4 Nanotubes are often found in bundles, 
ropes or fibers^ These three dimensional structures form 
through the van der Waals interactions between the tube 
surfaces, which also influence the properties of the bun- 
dles£i£ The orientation of tubes in a bundle can be com- 
pared to the stacking of graphene layers (e.g. bernal/AB 
stacked)) 10 ! 11 For nanotube bundles, however, we have 
to additionally consider the chirality of the tubes. The 
structural influence on the bundle properties is especially 
interesting for bundles of tubes that consist of only one 
chirality (monochiral bundles), as these have uniform 
properties needed for, e.g., electronic devices. 

The mechanical and electronic properties of bun- 
dles were already studied theoretically and experimen- 
tallyii2~— However, there is little knowledge about the 
influence of the structure (chirality) of the tubes on, e.g., 
the electronic structure of the bundle d 10 ! 18 Monochi- 
ral bundles, have not been experimentally produced yet, 
even so their production is expected in the near fu- 
ture^-^1 

Nanotube bundles and their electronic structure were 
studied with density functional theory within the local 
density approximation (LDA), which allows to consider 
the structure of the tubes in the bundle ^ 10 ' 18 The lo- 
cal density approximation, however, fails in modeling the 
van der Waals interaction between tubes. Another ap- 
proach is to use continuum approximations involving the 
Lennard- Jones potential. The continuum approximation, 
however, does not account for the specific configuration 



of the carbon atoms (structure) £2r2k Other approaches 
that consider the structure as well as the van der Waals 
interaction, do not reach the accuracy of density func- 
tional theory calculations jSLSI The van der Waals func- 
tional developed by Dion et al. was shown to be able to 
model bundles of carbon nanotubes in density functional 
theory 28 ! 29 At this level of theory studies considering the 
influence of nanotube chirality in bundles remain miss- 
ing. 

In this paper we study the properties of bundles con- 
sisting of carbon nanotubes of the same chirality and 
handedness by a van der Waals density functional. We 
calculated the binding energy, intertube distance, and 
electronic structure in dependence of the orientation of 
the tubes in the bundle. We find a particularly strong 
dependence on the orientation of the tubes in the bundle 
if the tubes are achiral and the chirality of the tubes and 
the bundle share the trigonal symmetry, meaning that 
the tubes have a Cg-axis. Achiral tubes with Cg-axis 
are metallic, e.g., (6,6) and (12,0) tubes. Pseudogaps 
open at the Fermi level in dependence of the orienta- 
tion of the tubes in the bundle. We find metallic be- 
havior for the AA-stacked configuration (0°) for Cg-axis 
armchair bundles, as well as for intermediate configura- 
tions which have glide reflection planes. The bundles 
of Cg-axis chiralities experience rotation barriers for the 
tubes that are induced by the configuration of the atoms 
on the tube surfaces. The barriers can be as high as 
AE = 10.5 meV/atom. Rotation barriers for tubes with- 
out the trigonal symmetry of the bundle are less than 
AE = 0.3 meV/atom. The properties of bundles made 
from non Cg-axis tubes have a much weaker dependence 
on the orientation of the tubes inside of the bundle, e.g. 
(8, 2) and (14, 0) tubes. 

This paper is organized as follows. We first describe 
the computational methods used for our study, Sec. II. 



We then discuss the van der Waals energies and intertube 
distances of bundles in orientational dependence in Sec. 
III. In Sec. IV we discuss the electronic band structure 
at the Fermi level for a (6,6) bundle as a function of 
tube orientation. Section V summarizes this work and 
discusses the conclusions in a general context. 



II. COMPUTATIONAL METHODS 

We used the ab-initio package SIESTA to perform our 
density functional theory calculations i^2r— To optimize 
the geometry of the isolated tubes we used the general- 
ized gradient approximation parameterized by Perdew, 
Burke and Ernzerhof, which is good in modeling the 
interactions within the tube.— We calculated the bun- 
dle properties (binding strength, intertube distance and 
bundle electronic structure) with the van der Waals den- 
sity functional parameterized by Dion et ai, which is 
especially good in modeling the van der Waals interac- 
tion between the tubes in the bundle.— Both calcula- 
tions used norm conserving nonlocal pseudopotentials. 34 
We balanced the computational time and accuracy, by 
choosing localized pseudoatomic orbitals with a doublc- 
C (DZ) basis set to describe the valence electrons. We 
chose cutoff radii for the s and p orbital of the carbon 
atoms with r s = 5.949 Bohr and r p = 7.450 Bohr. The 
mesh cutoff for the real-space integration corresponded 
to about 350 Ry. We used between 10 and 14 fc-points 
in the Monkhorst-Pack scheme 3 -^ to calculate the total 
energies in the /^-direction. The z-axis was chosen as 
the tube-axis. We performed all total energy calcula- 
tions with only one fc point in the x and y direction. The 
fc point sampling for the band structure calculations of 
the bundled tubes were 20 fc-points in x and y direction 
leading to a fc point sampling of 20x20x300 fc-points for 
(n, 0) zigzag tubes, 20x20x300 for (n, n) armchair tubes, 
and 20x20x200 for the chiral (8, 2) tube. In this paper we 
will only consider the band structure of the monochiral 
(6, 6) armchair tube bundle in detail. 

To obtain decent bundle structures we first optimized 
the lattice constant of isolated tubes by minimization of 
the total energy in dependence of the lattice constant. 
We then performed a geometry optimization for isolated 
carbon nanotubes within the generalized gradient ap- 
proximation to derive an initial tube geometry for our 
bundled tube calculations. We used the conjugate gradi- 
ent method to optimize the atomic coordinates to a max- 
imal force tolerance of 0.04 eV/A, which is adequate for 
our accuracy. The optimized isolated tube coordinates 
were then used to perform the van der Waals density 
functional calculations. 

For the van der Waals density functional calculations 
we placed a tube in a hexagonal unit cell with periodic 
boundaries, see Fig. [1] for the example of a bundle of 
(6, 6)-tubes. We calculate a bundle by an infinite number 
of tubes (bulk-bundle), which serves as a model of an 
inner tube / in a real carbon nanotube bundle. 
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FIG. 1: Ball and stick sketch of a bundle of seven (6, 6)-tubes 
in their high-symmetry orientation (0°). The hexagon in the 
middle (inner tube I) represents the unit cell for our bun- 
dle calculations. Lines connecting the middle points of the 
hexagons were added for clarity. The top carbon atoms are 
highlighted in blue, bottom carbon atoms are gray. The sym- 
metry for seven tubes is trigonal. The inner tube I interacts 
with six neighbors with intertube distance D. 



The total energy of the bulk bundle was minimized in 
dependence of the intertube distance to obtain the op- 
timal distance between the tubes of the bundle as pre- 
sented in Fig. [2l A Lennard- Jones potential is plotted 
with e and r m i n parameters fitted to our data. The 
Lennard-Jones potential shows a smaller width in the at- 
tractive region compared to the calculated data points. 
We did not perform a structure relaxation with the van 
der Waals density functional as this leads to energetically 
lowest states, which does not allow to consider energet- 
ically unstable orientations of the tubes of the bundle. 
Therefore, stress remains in the bundle systems depend- 
ing on the chirality and to a smaller extent on the orienta- 
tion of the tubes. The maximal forces in the systems are 
F m ax ~ 0.3 eV/A for armchair bundles, F max sa 1 eV/A 
for zigzag bundles, and F max w 3 eV/A for the (8, 2) and 
(12,6) bundle. Previous studies of tube-tube interaction 
found the distortion caused by the bundling to have lit- 
tle or no effect on the properties for tubes with diameters 
below 15 A£2£ 

The total van der Waals energy is defined as the dif- 
ference between minimal total energy of the bundle and 
total energy of the isolated tube (both calculated with 
the functional of Dion et a/.—), which corresponds to in- 
finite intertube distance. Dividing through the number 
of atoms of the unit cell yields the van der Waals energy 
per atom. We derive the intertube distance from the op- 
timized xy-unit cell length by subtraction of the diameter 
of the tube. The band structure was plotted for certain 
high symmetry directions. We used the whole fc-sampled 
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FIG. 2: Total energy in dependence of the intertube distance 
calculated for a bundle of (6, 6) tubes. The energy is normal- 
ized to the energy of an isolated tube, which corresponds to 
infinitely distant neigboring tubes. The data points were ob- 
tained by SIESTA. A Lennard- Jones potential is plotted with 
e and r m in parameters fitted to our data. 



zone to derive the density of states at an electronic tem- 
perature of 20 K. 



III. VAN DER WAALS ENERGIES AND 
INTERTUBE DISTANCES 

We start this section with a discussion of the symme- 
try properties of nanotube bundles. This discussion is 
followed by a study of the influence of the orientation 
of tubes in the bundle on the van der Waals energies be- 
tween the tubes of the bundle, as well as the energetically 
optimized intertube distances. At the end of the section 
we draw a conclusion from our results for the properties 
of mixed chirality bundles and monochiral bundles. 

Bundles were experimentally observed in a triangular 
lattice that contained up to hundreds of tubesji^i In our 
calculations we model the triangular lattice of the bundle 
by a hexagonal unit cell, compare Fig. [TJ A symmetry 
breaking occurs, if the individual tube and the bundle do 
not share all symmetry operations. The bundle structure 
has in general a D^h symmetry, this means it has symme- 
try axes (2 Ce, 2 C3, C2, 3 C'2, 3 C'2), mirror planes (3 
a v , 3 ad, crh), one inversion center (i), rotation- reflection 
axes (2 S3, 2 Se) and one identity element (E). The tubes 
have symmetry operations depending on their chirality. 
The level of symmetry breaking depends on the number 
of shared symmetry elements between the bundle struc- 
ture and the chirality of the individual tube. A tube in 
the bundle has to have the same atom configuration every 
60 on the circumference of the tube to share symmetry 
operations with the bundle. If the chirality and bundle 
structure share the symmetry, we say that the chirality 



TABLE I: Comparison of the intertube distance D, van 
der Waals energy per atom |_E^^|, and their variations 
A (_D, -E^dw^) f° r various chiralities (n,m), chiral angles 6 
and tube diameters d. The table presents the dependence 
of the fundamental properties on the orientation (rotation) of 
the tubes in the bundle with minimal and maximal values. 
Errors for intertube distances are ±0.02 A and for binding 
energies the errors are ±0.2 meV/atom. 
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has a C6-axis and symmetry breaking is lifted in high 
symmetry configurations (e.g. 0° in Fig. [1]). The high 
symmetry configuration has all symmetry operations of 
the Deh symmetry, especially mirror planes. 

We study the orientational influence on the properties 
of the bundles by simultaneous rotation of all tubes of 
the bundle in steps of 1 or 5 starting from to 60 °. 
We consider the (6,6), (12,12), (9,0), (12,0) and (12,6) 
Cg-axis chiralities. The (12,0)-tube structure, for exam- 
ple, has the same atomic configuration every 15 around 
the circumference (considering screw operations), as the 
full circumference contains 360 and there are 24 atom 
positions on the circumference of the tube. This also 
means that the (12,0)-tube has the same atomic config- 
uration at the angles 0°/360°, 60°, 120°, 180°, 240° 
and 300°, therefore it has a C6-axis. To study the ef- 
fect of intertube orientation on non Cg-axis chiralities 
we further study the chiral (8, 2)-bundle, the two zigzag 
bundles (8,0) and (14,0) and the two armchair bundles 
(5,5) and (10, 10), which have nearly no common sym- 
metry elements with the bundle structure. 

We present the minimal values, maximal values and 
variation in dependence of the orientation of the tubes 
(rotation angle) for the intertube distance D and van 
der Waals energy per atom for various chiralities 

and tube diameters, see Tab. U We only see small vari- 
ations for the values of non Cg-axis chiralities, e.g. for 
the bundle of (8, 2)-tubes. The intertube distance of the 
(8, 2)-bundle varies by AD w 0.03 A and Ef$? varies 
by s» 0.3 meV, meaning by 1% or less. Our results cor- 
respond well to results of the (10, 10)-bundlc, with our 
maximal difference between lowest and highest van der 
Waals energy per atom AE v dw = 0.3 meV and the ac- 
tivation barrier for rotations of AE vc iw = 0.15 meV re- 
ported previously*^ The activation barrier results from 
the reduction in symmetry from T)2h to C2h due to rota- 
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tionJ^ The properties of bundles made from non C6-axis 
chiralities show hardly any dependence on orientation. 

Comparing the intertube distance for the (8, 0) and 
the (14,0) tube, we find a variation of less than 1-2%, 
whereas the tube diameter increases by ~ 73 %. The 
van der Waals energy per atom decreases about 23 % for 
the same diameter comparison. The tube diameter has 
a strong influence on the binding strength, but does not 
influence the intertube distance. 

For the C6-axis chiralities, we receive a radically dif- 
ferent result for the dependence of the properties of the 
bundles on the rotation angle. We observe a variation 
of the bundle properties in dependence of the orienta- 
tion for all Cg-axis chiralities, compare Tab. [I] Starting 
with the zigzag chirality (9,0); The intertube distance 
varies by AD s» 0.18 A, corresponding to about 6 % 
variation, which is six times higher than the variation 
observed for non C6-axis chiralities and at least three 
times higher than the variation associated with tube di- 
ameter. Furthermore AE^yy ss 5.2 meV, corresponding 
to about 15 % variation, which is 15 times higher than 
the variation observed for non Cg-axis chiralities and only 
slightly smaller than the variation accounted to the tube 
diameter. For the (12,0) chirality the orientational de- 
pendence on the properties weakens, with AD « 0.15 A 
(w 4%) and AEf°™ « 3 meV (« 11%). The rotation en- 
ergy barrier AE^ff « 3.0 meV corresponds well to the 
barrier of 3.0 meV previously reported by LDA calcula- 
tions^ The influence of the orientation on the bundle 
properties decreases with increasing tube diameter limit- 
ing the occurrence of the special properties of the Cg-axis 
bundles to small diameter nanotubes. We suppose that 
the increase in interaction area lowers the influence of the 
local symmetry. The interaction area increases through 
the reduced curvature of larger-diameter tubes. 

The (6, 6) armchair bundle shows the strongest ori- 
entation dependence with AD ss 0.29 A (« 9%). The 
rotation energy barrier AEfjffi « 10.5 meV (« 32%) 
is larger than * 5.0 meV previously reported by LDA 
calculations. 10 Interestingly, the variations of the bundle 
parameters in dependence of the orientation are smaller 
for the zigzag bundles than for the armchair bundles of 
comparable tube diameter. 

For C6-axis bundles of chiral tubes (e.g. (12, 6)) we see 
a dependence on the van der Waals energy (AE^^ 1 « 
0.9 meV (~ 4%)) and the intertube distance (AD « 
0.04 A (w 1%)), however, the dependence is low. In 
chiral Cg-axis bundles alignment of neighboring tubes is 
impossible due to the handedness of chiral tubes. Neigh- 
boring tubes with same handedness have opposing sur- 
faces with different handedness leading to mismatched 
surface structures of neighboring surface atom layers. It 
is not possible to interchange handedness of neighboring 
tubes to generate alignment for all neighbors as neigh- 
bors with same handedness remain caused through the 
trigonal structure of the bundle. 

The intertube distance of achiral Cg-axis nanotubes 
depends on their orientation and varies by as much as 
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FIG. 3: Total energy for various intertube distances as a func- 
tion of rotation angle for a bundle of (6, 6)-tubes. The total 
energy depends on the intertube distance as well as the orien- 
tation of the tubes. The inset shows a side view sketch of the 
atomic configuration for the three most interesting rotation 
angles. The front tube has blue atoms, the tube in the back 
is in gray. 



10%. The binding strength of Cg-axis bundles is influ- 
enced through the variation of the intertube distance but 
at the same time by the orientation; see Fig. [3] The in- 
fluence of the intertube distance on the total energy has a 
strong impact for 0° (AA-stacked) and around 10° (AB- 
stacked), but is weaker around 4°. The orientational in- 
fluence on the binding strength of Cg-axis chiralities with 
small diameters (< 20 A) is of the same magnitude as the 
influence of the tube diameter on the binding strength in 
the bundle. 

We use the example of the (6, 6)-bundle to discuss the 
physical properties of a Cg-axis bundle. In this bundle 
the properties depend most strongly on rotational orien- 
tation. The van der Waals energy per atom of the (6, 6)- 
bundle as a function of rotation angle shows symmetry 
breaking behavior, see Fig. 01 Symmetry is initially Dg/j 
(0°) and reduces to Cg/j (loss of mirror planes) due to 
rotation. The high symmetry configuration of , corre- 
sponding to AA stacking in graphite^ -, is not energeti- 
cally stable, see Fig. [3] It has one of the highest intertube 
distances of 3.43 A. The magnitude of the van der Waals 
energy increases up to a maximum at ~ 10 °, correspond- 
ing to AB stacking in graphite^ -, see Fig. @] This is in 
contrast to previous LDA findings, where a configuration 
about 2.5° off the AB stacking lead to a maximum at 
about 7.5 ji£ The configuration at 10° has one of the 
smallest intertube distances observed in our calculations 
with 3.15 A. At 15 ° the rotation leads to an interesting 
configuration that has glide reflection planes, with a lo- 
cal binding energy minimum and a moderate intertube 
distance of 3.25 A. Further rotations only reproduce the 
behavior that is contained in the first 15 ° rotation, see 
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FIG. 4: Van der Waals energy per atom Evdw as a function 
of rotation angle for a bundle of (6, 6)-tubes. The tubes of the 
bundle are rotated starting from the high symmetry position 
(0 °, see Fig. [1]). The lines are cubic spline fits and show sym- 
metry breaking behavior. Global extrema occur every 10°. 
Minimum values occur at 10 ° + 30 ° • n and 20 ° + 30 • n with 
integer n. Global maxima occur every 30 ° starting from ° 
and local maxima occur every 30° starting from 15°. The 
insets show the two exemplary configurations of 10 ° and 15 ° 
rotation of the tubes in the bundle. The blue spheres mark 
the position of the top carbon atoms, the gray spheres mark 
the position of the bottom carbon atoms. 



Fig.m 

The error in our calculation can be estimated by com- 
paring two identical configurations, e.g. 10°, 20°, 40°, 
and 50 ° in Fig.UJ We obtain an error of ±0.2 meV/atom 
for the total energy. An error of ±0.02 A for the intertube 
distance can be estimated from the data in Fig. [3] 

The smaller activation barrier for rotations between 
15° and 20° has a value of AE = 1.8 meV/atom and 
the second, larger activation barrier for rotations be- 
tween all orientations (e.g. between 0° and 10°) is 
AE = 10.5 meV/atom. Rotations of a solid made 
from Cqq fullerenes were experimentally and theoretically 
found to be hindered below T « 260 K, implicating that 
a hinderance of rotations for nanotubes in a bundle might 
also be possiblei 14 i 38 ~ — 

At the end of this section we want to discuss our 
findings and the implications for mixed chirality carbon 
nanotube bundles. Most carbon nanotube bundles con- 
tain nanotubes of various chiralities (n, m). The binding 
strength and intertube distance depend in general on the 
orientation of the tubes in the bundle. But the orien- 
tational dependence of the aforementioned properties is 
suppressed by symmetry breaking induced by mismatch 
of the bundle and tube symmetry. We can conclude from 
our calculations, that the properties of bundles of mixed 
chiralities have a negligible dependence on tube orien- 
tation. This results from the non C6-axis chirality ob- 
servations, which show only a weak dependence on the 
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FIG. 5: Band structure along high symmetry directions of a 
bundle of (6, 6)-tubes rotated starting from its high symmetry 
configuration (0 °) in rainbow colors (color online). The band 
structure diagrams are normalized to the Fermi levels. 



orientation for the bundle properties (about or less than 
1% variation in orientation dependence). The tube di- 
ameter (curvature) otherwise has a strong influence on 
the binding strength, but does not influence the intertube 
distance. Recently Crochet et al. succeeded in producing 
nearly mono chirality single walled carbon nanotube bun- 
dles, highlighting the possibility to produce monochiral 
samples in the near future which would allow to experi- 
mentally test the theoretically predicted properties *3r— 
Only the monochiral C6-axis tubes preserve the symme- 
try in the higher symmetry configurations. Chiralities 
with C6-axis can exhibit special properties, e.g. binding 
energy and intertube distance are influenced in orienta- 
tional dependence in contrast to the tube diameter which 
only influences the binding energy. 



IV. ELECTRONIC BAND STRUCTURE AT 
THE FERMI LEVEL 



In this section we study the band structure, density 
of states, and pseudogap and their variation in bundles 
of metallic tubes belonging to the Cg-axis chirality at 
the example of (6, 6)-tubes. At the end of the section we 
compare our results to other theoretical and experimental 
studies and discuss the influence of rotation barriers on 
the electronic properties of bundles. 

Band structure diagrams of a bundle of (6, 6)-tubes 
for orientations of 0° to 15° show that the bands change 
their slope with orientation, see Fig. [5j The valence band 
maximum and conduction band minimum shift in energy- 
and fc-direction through the change in orientation. This 
is accompanied by a shift of the Fermi level, as can be 
seen between A and Y (all Fermi levels were normalized 
to eV). The Fermi levels shift through the changed 
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FIG. 6: Band structure along high symmetry directions of a 
bundle of (6, 6)-tubes starting from the position of the valence 
and conduction band crossing in rainbow colors (color online). 
The band structure diagrams are normalized to the Fermi 
levels. The valence band flattens by change in orientation 
while the influence on the conduction band is smaller and it 
mostly keeps its shape. The bands for the orientations with 
metallic behavior (without pseudogap) are highlighted. 



electron density of each orientation, which result from a 
volume change induced by the orientational specific in- 
tertube distances. For the high symmetry configuration 
(0°, Dg/i) the valence band and conduction band cross at 
about 43 meV above the Fermi level, which corresponds 
well to the value of about 70 meV calculated within LDA 
calculations; the Fermi level is shifted compared to the 
Fermi level of the isolated tube<2 The k z value for the po- 
sition of the valence band maximum and conduction band 
minimum is between k z = 0.60 • - for and k z — 0.66 • - 
for 15°, where a — 2.491 A is the lattice constant along 
the tube. Our fc z -values correspond well to the value in 
isolated tubes of k z = 2/3 • ^. There is no band splitting 
in the orientations at and 15°. The different par- 
ity of the bands allows the crossing of the bands. The 
band splits by rotating the tubes of the bundle as little 
as 1 out of the high symmetry position, which opens 
up a band gap of E g = 75 meV along the AT-direction. 
The mirror symmetry for rotated configurations is bro- 
ken, which leads to anticrossing. The band gap increases 
until « 7 and then becomes smaller again until it van- 
ishes at 15 °. The valence band flattens in dependence of 
the orientation for the high symmetry directions start- 
ing from the valence and conduction band crossing, see 
Fig. El 

The high symmetry configuration (0°) shows metallic 
behavior while an increase of rotation angle of the tubes 
in the bundle increases a pseudogap, which is maximal 
at w 7° , see density of states in Fig. [7] The pseudogap 
diminishes with further rotation until it closes at 15°, 
which can be accounted for by the symmetry of the con- 



figuration. The 15 configuration looses mirror planes, 
but has glide reflection planes, which keeps the configu- 
ration metallic, see Figs. El and [Jj Furthermore, the den- 
sity of states minima shift around the Fermi level from 
meV at the global minimum to 50 meV in dependence 
of the orientation, see right panel in Fig. [7] The shift of 
the density of states maxima leads to increased density 
of states at the Fermi level for orientations between 12 
and 14 , compared to the metal like behavior of the 
orientation. This can partly be understood by compari- 
son to the behavior observed by the band structure, see 
Fig. El 

The orientation has a smaller influence on the valence 
side than on the conduction side of the band structure, 
see Fig. [71 The first peak of the conduction side shifts 
between 1.07 eV and 1.22 eV through rotation; for certain 
rotation angles (e.g. 15°) it becomes a shoulder. The 
valence band shows only small variation for the first two 
peaks. 

We conclude that the most exciting properties of C6- 
axis bundles can be found in the band structure and den- 
sity of states. Achiral C6-axis bundles can be metallic. 
Certain orientations possess higher density of states at 
the Fermi level than the metallic configurations. The 
minima of the density of states meanwhile shift around 
the Fermi level in dependence of the orientation. The 
band structure shows a dive of the conduction band min- 
imum below the Fermi level along the high symmetry 
direction accompanied by a pseudogap opening which 
closes again for higher symmetry configurations that have 
mirror planes or glide reflection planes. 

We now want to compare our results to other theoret- 
ical studies. The density of states of the high symmetry 
and rotated (10, 10) monochiral nanotube bundles were 
shown to have very similar behavior, which agrees with 
our results i 14 ' 15 For C6-axis bundles, however, we find 
a strong coupling between the electronic states near the 
Fermi level and the rotational motion, which was sug- 
gested to lead to superconducting behavior previously.— 
In contrast to previous studies we find metallic behav- 
ior for achiral C6-axis bundles (e.g. (6,6), (12,12)) for 
multiple orientations^ We were able to show, that no 
energy gap, but a pseudogap opens in dependence of ori- 
entation of the tubes in the bundle in contrast to previ- 
ous studies, which showed an energy gap in the density 
of states for AB-stacked (6, 6) bundles as well as for the 
8° orientation^ Furthermore our results show a shift 
of the density of state extrema around the Fermi level, 
corresponding to previously reported results^ We also 
find, however, an increased density of states at the Fermi 
level compared to the metallic configurations for certain 
orientations. We find good agreement in comparison to 
experimentally derived density of states with the general 
trend of pseudogaps opening due to the bundling. 16 

At the end of the section we want to discuss the in- 
fluence of rotation barriers on the electronic properties 
of monochiral achiral Cg-axis bundles. At low temper- 
atures, certain configurations of monochiral Cg-axis chi- 
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FIG. 7: Density of states of a bundle of (6, 6)-tubes rotated starting from its high symmetry position (0°) in rainbow colors 
(color online). The density of states diagrams are normalized to the Fermi levels. The right panel shows the shift of the 
pseudogap minimum with respect to the Fermi level (0 eV), which can also be observed in Fig. [S] 



rality bundles can be stable, e.g. the orientation at « 10° 
for the (6, 6)-bundle. At room temperature thermal en- 
ergy is likely to lead to rotations and vibrations of the 
tubes in the bundle. A complex electronic behavior is ex- 
pected, as rotations (orientational changes) and absorp- 
tion and emission occur time averaged over the experi- 
ment. Especially the slope flipping behavior and shifts of 
the bands lead to a further broadening of the width of the 
bands, which further broadens the density of states and 
therefore also the optical-absorption bands for the Cq- 
axis chiralities. The high symmetry configuration with 
metallic behavior is the energetically most unstable con- 
figuration and therefore semi-metallic behavior has to be 
expected for bundles of metallic tubes. The 15° config- 
uration is a local total energy maximum and therefore 
also expected to be an unstable configuration, even so it 
is much lower in energy than the high symmetry config- 
uration. For larger diameter tubes the rotation barriers 
flatten, allowing metallic behavior in Cg-axis bundles of 
large diameter armchair tubes. 

V. CONCLUSIONS 

In summary we presented van der Waals energies and 
intertube distances of various chiralities of carbon nan- 
otube bundles in dependence of the orientation of the 
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